Effective mass of one 4 He atom in liquid 3 He 
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» ( , A microscopic calculation of the effective mass of one 4 He impurity in homogeneous liquid 3 He 

at zero temperature is performed for an extended Jastrow-Slater wave function, including two- 
^5 ' an d three-body dynamical correlations and also backflow correlations between the 4 He atom and 

the particles in the medium. The effective mass at saturation density, 7714/7714 = 1.21, is in very 
good agreement with the recent experimental determination by Edwards et al. The three- particle 
correlations appear to give a small contribution to the effective mass and different approximations 
for the three-particle distribution function give almost identical results for 7714/7714. 
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■ Recently, Edwards et aM have used the Zharkov-Silin 
I Fermi liquid theory of dilute solutions of 4 He in normal 
. liquid 3 He to determine the chemical potential (/xi) and 
t^J- ' the effective mass (m|) for the limiting case of one 4 He 
OS . impurity. The experimental input data for this analy- 
' sis were the recent low temperature measurements of the 
phase separation by Nakamura et ala The same theory 
£3 ■ was previously applied to older experimental data by La- 
1 , I heurte and Saam,Hul and their predictions for /Z4 and m| 
differ notably from those of Ref. |l| At zero pressure, EcL 
wards et al. report /14 = —6.95 K and m\/rrii — 1.1, □ 
Zj , whereas /14 = —6.60 K and 7714/7714 = 4.5 in Rcfs. |]||. 
' The disagreement was attributed to the fact that the 
. ,-h , validity of the Fermi liquid theory is ensured at tempera- 
^ ■ tures below 0.1 while the results of Ref. |,| were ob- 
5h I tained from experimental data at T > 0.5 K. Although 
. . . 1 the difference between the two experimental values, of \i± 
is rather small, our recent microscopic calculationsO seem 
to support the more bounded result of Edwards et al. On 
the other hand, the experimental determinations of the 
impurity effective mass are appreciably different. At this 
point, it is clear that a fully microscopic calculation of 
m\ would be very enlighting. 

In the present work, we evaluate the excitation spec- 
trum and the effective mass of a 4 He impurity in liquid 
3 Hc using a trial wave function of the type 

tf„(k) = pfl(k)tf , (1) 

where ^0 is the ground-state wave function of the 3 He 
medium plus one 4 He atom and pbO<-) is an excitation 
operator defined as: 

p B (k) = Pl (k) F B , (2) 



1 



where p/(k) = exp(ik • ri) describes the impurity travel- 
ling through the medium as a plane wave of momentum 
k and the correlation operator 

A 

F B = ]Jf B (k,r Ii ) (3) 

i=l 

incorporates backflow correlations between the impurity 
/ ( 4 He atom) and the A 3 He atoms of the bulk. Back- 
flow correlations have proved to be relevant for a realistic 
stucUj af the effective mass of one 3 He impurity in liquid 
4 He.B~0 They play also an important ralajn the evalua- 
tion of the binding energy of pure 3 He.mL9 
The Hamiltonian of the system is written as 

H{A + 1) = H{A) + Hi{A + 1) , (4) 

where 

^) = -^E v '+X>M ( 5 ) 

i=l i<j 

is the Hamiltonian of the pure 3 He background, and 

2 A 

ffi(A + l) = -A-V? + ^y(r i/ ) (6) 
4 i—i 

are additional terms related to the impurity. 

The variational approach starts with the choice of. the 
trial wave function ^o- As in our previous paper p we 
take an extended Jastrow-Slater wave function for the 
A + 1 particles: 

y =FjF T cj>(l,...,A) , (7) 

where 0(1, A) is the Fermi gas wave function for the 
A 3 He atoms, Fj is a Jastrow correlation operator em- 
bodying two-body dynamical correlations 

^=n/ M yn/ M (^). (8) 

i<j i—1 

and the triplet correlation operator Ft is written as: 

A A 

F T = J[ e-^^IJe-^/ 2 , (9) 

i<j<k i<j 

with 

q a jk = ^£,(r aj )cl(r ak )r a: j ■ r ak . (10) 

eye 

Here, ^2 cyc denotes a summation on the cyclic permu- 
tations of the indices a,j, k, where the index a may ei- 
ther represent a 3 He atom or the 4 He impurity. The 
triplet correlations have been found to be very impor- 
tant to properly describe the equation of state of the 
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pure phased and to calculate the chemical potential of 
the 4 He impurityE 

The expectation value of the Hamiltonian with respect 
to ^ v (k) is given by 



*0 



(*o I *o> 



(11) 



where we have taken advantage of the unitary character 
of the excitation operator ps(k). If one ignores backflow 
correlations between the impurity and the medium, by 
assuming /9b (k) = p/(k), a simple parabolic spectrum 
for E v {k) is obtained: 



E v (k) = E% 



2 U2 



h z k 

IrriA 



(12) 



where Eq = {^q\H In this case, the 4 He effective 
mass is equal to the bare mass, pointing to an excessively 
simple choice for the excitation operator. 

A better ansatz is given by W„(k) of Eq. which 
explicitly contains backflow correlations between the im- 
purity and the 3 He atoms. In particular, the Jpackflow 
correlation operator has been taken of the fornxra 



/ /; = \\ cxp [ik ■ vn r]{rn)] 



(13) 



As the interatomic potential depends only on the rel- 
ative distance between the atoms, it commutes with 
ps(k). Therefore, in Eq.(|ll]) it is only necessary to con- 
sider the commutator with the kinetic energy operator. 
After some integration by parts, a generic contribution of 
the kinetic energy operator to Eq.(|ll|) may be expressed 
through the following relation: 



X a — 



Pb NI,Pb} 



with 



X p = - 



(*o I *o) 

(V a P B ) (V Q ps)|*o) 



xi 



and 



XI = 



(*o I *o> 



(14) 



(15) 



(*o I *0 



V Q p\ (V a p B ) (FjF T ) 



(*o I *o> 



(16) 



where the subscript a labels the generic particle. The 
arrows indicate in which direction ( right or left) the 
derivatives are acting. 

The first term X? ( |l5| ) is analogous to the expres- 
sion obtained in the case of the 3 He impurity in liquid 
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4 He.B! The second one X% (|l6|) directly originates from 
the Fermi character of the 3 He medium, as it comes from 
the kinetic energy operator acting on the Slater deter- 
minant 4>. Clearly, a — I docs not contribute to X%. 
Moreover, by inspecting the cluster expansion of X£, it 
results to be strictly zero. In fact, for direct cluster terms, 
where the a-particle is not exchanged, each of the two 
pieces of X% is zero, after summing over the momentum 
carried by a. Terms in which a is exchanged cancel be- 
cause each of them gives the same contribution in both 
pieces. 

The explicit expression for the impurity single-particle 
excitation energy, measured with respect to the 4 He 
chemical potential, is then 



e k = E v {k)~E% = 



2 1,2 r 



hrk 

2?7l4 



1 + e 2 + 



m,4 



e 3 



. (17) 



where 



('2 



dv ij 9?j 



2 Vij 



(18) 



j (2) 



(19) 



and 



e 3 = P 2 J dvijdrik gf) k 



VljVlk 



(20) 



fi is the reduced mass (fi 1 = 



_1 ) and g*j) and 



(3) 

9i/k are ^ ne two- and three-body distribution functions 
between the impurity and the 3 He atoms of the medium. 
They are the only quantities carrying information about 
the antisymmetry of the 3 He bulk. It is worthwile to 
remind that, by changing with 7713 in Eq. ([l7|), one 
recovers the expression for the reverse problem of one 
3 He impurity in liquid 4 He (Eq. 2.26 of Ref. § ), with 
the obvious substitution of the appropriate distribution 
functions. 

All the calculations presented in this paper have been 
performed in the framework of the so called Average Cor- 
relation Approximation (AC A). In this approximation 
one considers the same dynamical correlation functions 
for all the pairs and triplets in the system, not distin- 
guishing between the two isotopes. This assumption re- 
lies on the fact that the interatomic potential is the same 
for all the pairs. The drawbacks of the AC A in the evalu- 
ation of the chemical potential of the 4 He impurity have 
been extensively discussed in Ref. ||. _ 

We have used the interatomic Aziz potentiate and the 
two-body correlation factor f(r) has been taken of the 
McMillan type: 
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f(r) = exp 



1 (i 

2 \r 



(21) 



The variational parameter b has been fixed by means 
of a numerical minimization of the energy of pure liq- 
uid 3 He. The value b = 1.15 a (a = 2.556 A), de- 
termined at the 3 He experimental equilibrium density 
(P(T P = 0.277 er~ 3 ), has been used for all the densi- 
ties. The function £(r) of the triplet correlation Jl0| ) 
has the same— parametrized form used in pure phase 
calculationsoilil: 



A* exp 



tot 



(22) 



The density dependence of the triplet variational param- 
eters is neglected and the optimum values at p^ p have 
been used everywhere. These_yalues are A t = —0.75 ct~ 2 , 
r t = 0.85 g and w t = 0.45 cr.B 

The distribution functions have been computed by us- 
ing the Fermi HyperNetted Chain (FHNC1 technique, in 
the so called FHNC/S(T) approximation!] to take into 
account the elementary diagram ( and triplet) contri- 
butions. As reported in Ref. |5| the chemical poten- 
tial provided by the variational wave function a t 
the FHNC/ST saturation density (p = 0.252 cr" 3 ), is 
-6.60 K. 

The function rj{r) (jlffi— adopted for the backflow cor- 
relation, is of the forrrlroO 



r)(r) = A exp 



w 



(23) 



In our case, the backflow parameters Aq — 0.2, ro = 0.8 a 
and u>o — 0.375 a, taken from Ref. |l^, are used at all the 
densities. 

As the single-particle spectrum ( |T7j ) is quadratic in k, 
the effective mass is given by 



mA [/51 = 

m iJ 1 + e 2 + (m A /n) e m + , 



1 



(24) 



where (3 labels the approximation used in evaluating the 
three-body distribution function. 

Table I reports the effective mass obtained at the ex- 
perimental equilibrium density p e $ p in different approxi- 
mations. Also given is the value of the two-body contri- 
bution {m%/m^ 2 \ i.e. taking ef = in Eq. (|J). 

The Jastrow (FHNC/S) and the Jastrow plus Triplet 
correlation (FHNC/ST) models give nearly the same re- 
sults. The three-body distribution function has been 
evaluated in the Kirkwood superposition approximation 
(KSA), in the convolution approximation (CA) and in- 
cluding the Abe terms (KSA+ABE)£i As shown in the 

Table, the three different approximations to g\- k give 
very close results. 

The effective mass at p e xp in FHNC/ST, with the in- 
clusion of the Abe diagrams, turns out to be 1.21. This 
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result is in very good agreement with the most recent ex=. 
perimental determination m\jm^ = (1.1 +jQi4 / — 0.1).u 
As it has been pointed out by Leggettji-a the effec- 
tive mass of one impurity in a Fermi liquid is always 
larger than the bare mass. In FHNC/ST, e 2 = -0.37, 



e m = 0.08 and the three-body term (|20j), in all approx- 



imations, is very small, ~ 0.01. Therefore, the de- 
nominator of Eq. ( p4| ) is smaller than unity, providing 
an effective mass larger than one. 

The HNC/S results have been obtained by setting 
<f>o = 1 in the wave function, i.e. by treating the 3 He 
as a bosonic fluid. The comparison with the FHNC/S 
results indicates that, at this density, the influence of the 
Fermi character of the medium on the calculation of the 
effective mass of the impurity is nearly negligible. 

The density dependence of the calculated effective 
mass is reported in Table II and it is also shown in Fig. 
1 (full triangles). As one can see, the effective mass in- 
creases linearly with density. 

It is also interesting to compare our results for m\jm^ 
with the effective mass of one 3 He impurity in liquid 4 Hc. 
As there are not exchange or spin correlations between 
the two isotopes, one expects the effective mass to be 
driven mainly by the density. To deeper explore this 
hypothesis we have plotted in Fig. 1 the density depen- 
dence of the effective mass of a 3 He impurity in-Jiquid 
4 He, for both the experimental data (full circles)Ej and 
the theoretical estimates, .obtained by using backflow cor- 
relations (empty circles). a The density dependence is in 
both cases approximately linear. Although the slopes are 
different, the extrapolated values of m^/m^, at the 3 He 
saturation density, are similar and close to the present 
evaluation of m^/rn^ As it has been mentioned before, 
by taking the proper mass factor 7713//Z in front of e m , 
Eq. ( |l7|) is approximately valid for one 3 He impurity in 
liquid^3e, since the differences_in the distribution func- 
tions in the two cases are small. E£l In fact, if one performs 
a calculation at p^ p using the mass factor 7713//U, then 
the HNC/ST result (shown by an open diamond in Fig. 
1) coincides with the backflow extrapolated value. The 
small difference between this value and the proper result 
(full triangle) of m\/m± arises almost completely from 
the different mass factors in front of e m . The difference 
practically coincides with the estimate 
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mi J 1 + e 2 + (m 3 /^) e m + e 3 

= 0.073, (25) 

1 + e 2 + (m 4 //i) e m + e 3 

obtained by considering the same distribution functions 
in the two systems. Assuming a linear extrapolation, the 
experimental value for m^/m^ (solid line) is close to the 
backflow extrapolation (long-dashed line) at p e ^ v , point- 
ing out that possible perturbative corrections, beyond the 
backflow terms, are small at this low density. These cor- 
rections have been evaluated! in Correlated Basis Func- 
tion theory (CBF) for 7713/7713 in 4 He. They result to be 
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about 0.5 at the 4 He saturation density (p = 0.365 <J~ 3 ) 
and rapidly decreasing with the density. 

It is worthwile to notice that the effective mass of a-iHe 



much larger than the effective mass of the 4 He impurity. 
The statistics and the spin effects, which are suppressed 
in the case of the 4 He impurity, appear to be the main 
responsibles for this difference. 

To briefly summarize, we have calculated the effective 
mass of one 4 He impurity in liquid 3 He by using back- 
flow correlations. These correlations provide for an ac- 
curate description of the He impurity spectrum at low 
momenta. Our results support the new experimental de- 
termination of Edwards et alB and.-are far from the pre- 
vious result of Laheurte and Saam.ua 
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FIG. 1. Density dependence of the impurity effective mass. 
Full and open circles are respectively the experimental data 
and the backflow results for 7713/7713. The open diamond is 
the backflow result for 7713/7713 at pg Xp , indicated by an arrow. 
The full triangles are the backflow results for mX/rrn. The 
experimental result is plotted as a full diamond with its error 
bar. The lines are linear fits to the corresponding points. 

TABLE I. 4 He effective mass at p£ v = 0.277 cr~ 3 in differ- 



ent approximations. 





HNC/S 


FHNC/S 


FHNC/ST 


(ma/m 4 ) 121 
(m 4 /m 4 )^ SA i 
(ml/ mi ) lCA] 
{ml/m^ KSA + AB ^ 


1.209 
1.197 
1.187 
1.197 


1.221 
1.208 
1.196 
1.206 


1.225 
1.213 
1.200 
1.210 


TABLE II. 4 He effective mass 
FHNC/ST approximation. 


as a function of density in 


P(*- 3 ) 


0.253 


0.277 1 


D.300 0.330 


(m 4 /m 4 ) 121 
{ml/m^ KSA+AB ^ 


1.20 
1.19 


1.22 
1.21 


1.25 1.28 
1.23 1.26 
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